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ABSTRACT 

The large N c limit of Yang-Mills gauge theory is the dynamics of a closed gluonic chain, but this 
fact does not obviate the inherently strong coupling nature of the dynamical problem. However, we 
suggest that a single link in such a chain might be reasonably described in the quasi-perturbative 
language of gluons and their interactions. To implement this idea, we use the MIT bag to model 
the physics of a nearest neighbor bond. 
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The sum of (planar) diagrams which survives 't Hooft's large N c limit of SU(N C ) gauge theory 
]|] can be interpreted as the quantum dynamics of a closed chain of gluons |2| . The limit establishes 
a cyclic ordering of the gluons in a color singlet state, and only nearest neighbors on each cycle 
interact. As long as one stays with the continuum theory, there are no further simplifications pos- 
sible in the asymptotically free situation: dimensional transmutation eliminates the dimensionless 
coupling, trading it for a scale, which sets the mass gap if there is one. In that case there are no 
free parameters in the theory, and its single particle spectrum is knowable only in the context of 
a complete solution; it is an all or nothing affair. The case of J\f = 4 supersymmetric Yang-Mills 
theory is an exception because that theory is superconformal and hence scale invariant. Then the 
't Hooft coupling N c g 2 is a free parameter, which can be varied in the hope of further simplifica- 
tion. This is the setting for the dramatic developments of the past year surrounding Maldacena's 
proposal that the large coupling limit is equivalent to supergravity on a five dimensional Anti de 
Sitter space-time || |3|, ||] . 

A quantitative understanding of the single particle mass spectrum at N c — ► oo in a non scale- 
invariant gauge theory, such as QCD, remains an apparently formidible problem, although there are 
suggestions and conjectures that some form of the Maldacena conjecture might yield a solution Q. 
Alternatively, one can introduce a cutoff of some type, thereby temporarily adding a parameter and 
providing the theory with a flexibility for further simplification, albeit at the cost of straying from 
the continuum limit. The existing proposals for the QCD/5 dimensional string theory connection 
a la Maldacena effectively do this [Q]. The large N c limit of lattice gauge theory [||, [| is a more 
conventional version of this approach. Yet another is our early work with a light-front (infinite 
momentum frame) description, in which the P + carried by each gluon is discretized [10, |ll|, In 
this latter work the large coupling limit in the presence of this cutoff leads to a string interpretation, 
although the physical nature of the resulting string is highly sensitive to the details of the cutoff 
procedure. 

In this paper we put aside attempts at a detailed quantitative solution of large N c gauge theory, 
and instead aim at an understanding of the qualitative features of the gluon chain dynamics which 
will play a crucial role in a more complete treatment. The overriding question, of course, is whether 
the gluon interactions are sufficiently attractive to form a bond between each nearest neighbor pair 
in the chain. After that there is the issue of the fate of the gluon spin. A potentially disastrous 
outcome would be a ferromagnetic tendency for the spins to align yielding low mass glueballs 
with enormous spin. The spin-spin interaction between neighboring gluons should have an anti- 
ferromagnetic character. 

With regard to the first point, it is certain that the perturbative, renormalization group im- 
proved, Coulomb force N c a s (R)/2R 2 ~ 3tt / (11R 2 In RA) as R — > is attractive between neighbors 
due to the color structure imposed by the large N c limit. However, to bond a pair of massless 
gluons, the coefficient of 1/R 2 must be of order 1 (in order to balance the gradient of the kinetic 
energy of order 1/R), which goes beyond the validity of perturbation theory. As R increases from 
small values, this coefficient grows, and it is plausible that at some critical R c the attraction is 
strong enough to bond, although this speculation is out of perturbative control. However, if bond- 
ing of a neighboring pair is established, a chain of many gluons can then be built from similar 
bonds forming between each nearest neighbor pair. An appealing feature of this scheme is that it 
does not require extrapolating the basic perturbative gluon dynamics of any given pair too much 
beyond R c where a s = 0(1). 

We are close here to the philosophy of the MIT bag model of low mass hadrons 13, 14, 15]. 
In principle, the complete description of a light single glueball state should involve all scales, and 
there would be no a priori reason to expect the bare gluon interactions derived from the classical 



1 



Lagrangian to have much relevance. However the phenomenological success of the MIT Bag model 
indicates that this gloomy expectation is too pessimistic. In that model the large distance effects on 
a given single hadron state are lumped into a boundary condition on the fields which are confined 
to a cavity (for simplicity, usually but not necessarily taken spherical) with a state dependent size 
R. The fields within the cavity are then treated perturbatively, and the color magnetic part of the 
interactions yields a qualitatively correct pattern of hyperfine splittings within a multiplet with a 
fixed number of partons; for instance, this model works very well for the 3 quark baryonic states. 
The actual measured size of the splittings unfortunately requires the coupling to be uncomfortably 
large: N c a s /7r ~ 2. Nonetheless the idea that confinement effects are flavor and spin independent 
and that the spin dependent forces are well described by basic tree level interactions seems a 
reasonable working hypothesis for our qualitative analysis. 

We therefore examine the magnetic interaction between two neighboring gluons on a glueball 
chain. Since it will be important to compare the size of the gluon exchange contribution to the 
quartic contribution, we confine the two gluons to a spherical cavity.Q This provides an infrared 
cutoff in space, and also takes rough account of the interaction of the two gluons with the environ- 
ment of the rest of the chain. The first step in this analysis is to find the modes of noninteracting 
gluons in the cavity. This is a standard textbook problem (see for example chapter 16 of ]lq]). One 
simply solves Maxwell's equations with time dependence e~ %wt . 

V x E = iuB V x B = -iuE V • E = V • B = 0. (1) 

Then 

(-V 2 -u 2 )E = (-V 2 -uj 2 )B = 0, (2) 

the solutions of which are, of course, ji(u>r)Yi m (9, <p). To solve the divergence conditions, it is 
convenient to divide the modes into (TE) and (TM) modes: 

E£? = Kji(ur)LY lm B TE = —VxE TE 

ILU 

B™ = Kji(ur)LYi m E™ = —V x B™ (3) 

where we have taken only the regular solutions at r = 0. 

The frequencies uj are determined by the color confining boundary conditions of a perfect mag- 
netic conductor: E norm = B tan g = 0. The lowest non-zero solution for uj occurs for the (TE) 
case at I = 1, with the boundary condition zoj' 1 (zo) = —ji(zo) the lowest solution of which is 
ujqR = zq ~ 2.7437. If the cavity were an electric conductor, the boundary conditions would have 
the roles of electric and magnetic fields reversed, but we stick to the color confining conditions in 
this article. Finally, we note that we may take the vector potential of each mode to be simply E/iuj 
in which case it automatically satisfies the Coulomb gauge condition V • A = 0. It is clear from the 
construction of the modes that the angular momentum carried by each mode is just I. 

As already mentioned the lowest frequency mode is the (TE) case with 1 = 1: 

K K I 3 

A- m = - — ji(uJ r)LY lm = - — ji(u Q r)\ —r x e m (4) 

IUJq IUJq V 47T 

* The elements of the following calculations were drawn from unpublished bag model calculations I did in the mid 
seventies for the lightest glueball states with spin and 2. Some of my results were cited in Jjfl[ |. Similar calculations 
were independently done in . We repeat some of that early work here because the new context requires adaptations 
which are most easily explained within a complete self-contained account. 
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where the polarization vectors are e\ = e*_ 1 = i,0)/y/2 and eo = (0, 0, —i), and K is determined 
by the normalization condition 

K 2 r R 



f d 3 xA* m -A m = 2—j r 2 drj 2 (. 
J oj Jo 



K 2 

= —R 3 (jj(z )+ji(z )-3j {z )j l {z )/z ). (5) 

From the above considerations we know that these modes carry angular momentum j = 1. We 
would like to put two gluons in this lowest mode and then study the spin-spin interaction induced 
by the QCD interactions. At zeroth order the total energy of the system is 2ujq = 2zq/R. This 
energy level is 9- fold degenerate with j = 2, 1,0. We are next interested in how the interactions 
resolve this degeneracy. 

To find out, we set up cavity perturbation theory, in which the zeroth order fields are precisely 
the modes we have just discussed. We therefore expand the vector potential in a complete set of 
normal modes: 

A(x,t) = ^-^[a n a n (x)e-^ + ai<(x)e-i (6) 

where oj n is the frequency for the nth mode, and the a n are just the vector potentials of the normal 
modes, normalized according to 

J d 3 xa n (x) ■ a m (x) = 5 nm . (7) 
They also satisfy the completeness relation 

]T a*(x)a£(y) = (hi - ^?)5(x - y), (8) 

n 

where the delta function includes a transverse projector as required for Coulomb gauge. The gluon 
propagator has the representation 



i(0\T[A k ( X f a My) s 7 ]\0) 



<*n( x )«n(y) + «n*( x )«i(y) 



a 7 i 2n e ^ 

X S X (3 [ du^-iu,(x°-y )X- «»(x)o4*(y) (q\ 
a ^J 2tt 2_ ie W 



where we made use of the fact that J2u>„ /ixed a n( x ) a n (y) can be shown to be real. 

We are interested in the spin-spin interaction which has its origin in the magnetic part of the 
Hamiltonian, 



Hmag = I CPx^TrF^ 



= J d 3 x^{Tr(d k Ai - diA k ) 2 + 2igTi(d k Ai - diA k )[A k , A t ] - g 2 Ti[A k , Ai} 2 } 

= J d 3 x{^Tr(d k Ai) 2 + igTvd k A t [A k , A t ] - ^Tr A k At[A k , Mh (™) 



3 



where in the last line we have committed to Coulomb gauge V • A = 0. The lowest order splittings 
are 0(g 2 ), but of course the cubic terms will only enter at this order in second order perturbation 
theory. We first evaluate the contribution of the quartic terms for which first order perturbation 
theory suffices: we simply must evaluate the matrix elements of the quartic terms on the zeroth 
order states and diagonalize the resulting matrix. 

We choose the color structure of the zeroth order states corresponding to a single link in a 
glueball chain: 

\mi,m 2 )i = -j={4 ni a\ n2 f a \G). (11) 

Note that this state is in the adjoint representation of SU(N C ), as is appropritate for a single link 
in a gluonic chain. Contrast this with a typical bag calculation which is restricted to color singlets 
only. The operators in this state each create a gluon in the lowest cavity mode. We have suppressed 
all labels except spin mi^ = ±1, 0, and we recall that they are matrices in color space, a, (3 in ( |TT| ) 
being color matrix indices. The norm of this state is unity in the limit N c — > oo by virtue of the 
factor 1/y/Nc. Next we evaluate the matrix element in the limit N c —* oo: 

_„2 2]y 

^(m' 1 ,m' 2 \^-TTA k A l [A k ,A l ]\m 1 ,m 2 f a = 

_i_ k* k* I I _ n k* l# I k ] /'19"\ 
a m' 1 a m' 2 a m 1 a m 2 z,a rn' 1 a m' 2 a mi a m2l \ LZ ) 

Here the are just the mode functions for the lowest cavity mode contribution to A k (see @, (||)). 
Inserting these explicit expressions and integrating over space, we obtain the quartic contribution 
to the energy shift matrix (stripped of color factors) 

(m'i , m' 2 1 ^quartic \m\ , rm) 



The two terms in square brackets are commuting 9x9 spin matrices. The first has eigenvalues 
— 1, 1,2, where a spin 2 quintet belongs to —1, a spin 1 triplet belongs to 1, and a spin singlet 
belongs to 2. The second matrix has respective eigenvalues 0,2,0, so it only serves to raise the 
triplet two units yielding the pattern —1,3,2 for respectively spin 2, 1, and 0. This inverted level 
structure means that the quartic interaction yields ferromagnetic spin-spin couplings. We next turn 
to the contribution of the cubic terms for which we must use second order perturbation theory. 

To do second order perturbation theory we expand the interaction picture evolution operator 
to second order 

Te -i^ dt'Hj{t') =I _ { F dt > Hl (t') - - f dt! [ T dt" T[H I (t')H I (t'% (14) 

Jo 2 Jo Jo 

and consider only the contribution of the cubic terms of Hj in the last term. To extract the 
second order level shift we take matrix elements between the zeroth order states and isolate the 
coefficient of — %T for T — > oo. Since we are stopping at order g 2 , terms with higher powers of T 
will not appear, so this procedure gives the energy shift unambiguously. We will also exploit the 
simplifications of the limit N c —>oo. Four of the six factors of A in the last term contract against 
operators in the states, leaving the last contraction to produce a factor of the gluon propagator 

D kl (x,y,x° -y°)6 5 jP = i{0\T[A k (x)^My) s J\0) (15) 



N c g 2 3K 4 

167T UJ, 







R 



r 2 drjf(u:or) 



[(em! x e*. ) • (e m2 x e*. ) + (e mi x e m2 ) ■ (e*, x e*, )] (13) 
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which satisfies: 

(di-V 2 )D kl = 5Ux-y). (16) 
The cubic term in H mag can be written 

H mag , uUc = ig Jd 3 xTrd k A l [A k ,A l ]. (17) 

Inserting this into the second order term in the evolution operator, and contracting out a pair of 
gluon fields in all possible ways leads to 

--/ dt' [ d^ /, T[ J ff / (^ , )i^/(^ ,, )] - / d*x f d 4 yD kl (x,y)Tr T^MMy)], (18) 

* Jo Jo * Jo Jo 

where the current operator J k (x) is found to be 

J k {x) = At^Ai + [A k ,V- A] + 2VA k • A - 2A • VA k . (19) 

A contraction which involves d k Ai requires an integration by parts, but then the surface term is 
proportional to which is proportional to [A norm , At an g] and so vanishes for both electric 

(Atang = 0) and magnetic (A norm = 0) conductor boundary conditions. 

When (|l^) is sandwiched between two states of the same energy, a factor of T arises as T —* oo 
because then the integrand is only a function of x° — y°, but x°,y are integrated independently 
over the large interval T. Thus for Ef = E{ we have for the coefficient of — iT as T — > oo, 

(f\AE\i) = £ J^dx ! d 3 xd 3 yD kl (x,y)\ y0=0 (f\Tr T[J fe (x)^(y)|,o =0 ]|i). (20) 

Taking \i), \ f) to be of the form (11), we find that in the large N c limit there are only two essentially 
different contraction patterns that survive. Note that the equivalence of k,x to l,y means that we 
can insist that the first creation operator in \i), say, contracts with one of the two fields in J(y) 
provided we multiply by 2. After that contraction, which can be taken in two distinct ways, the 
remaining contractions are all uniquely determined when N c —>oo. There is an annihilation graph 
in which both operators in \i) contract against the fields in J{y). And there is an exchange graph 
in which each current contracts with one operator in \i) and one in |/). 

We first evaluate the exchange graph. Since the two gluons in our chosen states have the same 
mode frequency, the currents coupling to the propagator are time independent so the x° integral 
can be immediately done: 

f dx°D kl {x,y)\ y0=0 =^^[a£(xR*(y) + a£* (x)a^(y)] = G«(x,y). (21) 

J n Zljj n 

Clearly G satisfies the static Green function equation 

-V 2 G H (x,y)=^(x-y). (22) 

The matrix element reduces to 

(f\AE\i) exchange = N c g 2 5^5^, J d?xd 3 yG kl {^y)J k 2 {^)jhy) 

= N c g 2 5«'5P, j ' d?xji(*)A k (*) (23) 
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Here J is the matrix element, stripped of color factors, of J between one gluon states: 
+ a mu iV ■ a^, - V • a mi Q; ;ii + 2Va mi)J • a^, - 2a mi • Va^, ;) 

^fe( X ) = ^(«^r4«m 2 ,r 

+ «mi,,fc v ■ "ma - V ■ a^a,^,* + 2Va^ jfc • a m2 - 2a^ • Va m2jfc ) (24) 

Notice that ,7 fc 2 is the same as —J^ with the substitution mi, mi — ► m2,m' 2 . 

We do not need a closed form expression for G^. Instead we can work with Ak(x) = / d 2 yGkiJi, 
which satisfies 



- V 2 A k = J* = -^^^flr x (e mi x <,,)]*, 



(25) 



together with the magnetic conductor boundary conditions A n0 rm = and (V x A)tang = 0. The 
r.h.s. doesn't require the transverse projector because V • J 1 = 0. This equation is easily solved 
by the ansatz A = A(r)r x (e mi x e* , ), whence A(r) satisfies: 



A n , 4 y _ 9if 2 ji^prf 
r 47rojg 2u;o^ 2 
This last equation can be directly integrated to obtain 

m 2 



A{r) 



87TCJq 



(26) 



(27) 



where the integration constant is fixed by the boundary condition RA'(R) + 2A(R) = 0, which 
assures the vanishing of the tangential magnetic field. In terms of these quantities the level shift 
becomes 



(f\AE\i) exchange = N c g 2 8Z'6i^- r 2 dr4(r)ji(a;or) 2 (e mi x e* m , ) • (e m2 x e*,) 

Wn Jo 1 2 



3i£ 2 

'o 

9if 4 /" R (ir 



ii(w r) 2 (e mi x e*,)- (e m2 x e* ,) 



rr r 3 i-R 

-2 J r' 2 dr'ji(uQr') 2 - r' 2 dr' 'ji(u r') 



(28) 



To simplify this integral we note that 

X (z)= \ Z z' 2 dz'h{z) 2 
Jo 



■[zji(z) 2 + zjo(z) 2 - 3jo(z)j 1 (z)}. 



Putting z = zq gives the integral required to determine the normalization constant 2{K 2 / ujq)x{zo) = 
ujq. Numerically x( z o) ~ 1-1385. Stripping off the color factors, we have the level shift matrix for 
the exchange graph: 



(m' 1 ,m' 2 \A exchange \m 1 ,m 2 } = 



N c g 2 9lo 



16vr x(z ) 2 

N c g 2 uJo , p 



T° dz ■ f \2 f \ 
- / — Ji(z) 2 x(z) - 

Jo z 

( e mi X e 



4*g 



-0.4958- 



lg7r 5 a V ( e mi x e^,) • (e m2 x e^) 



(29) 
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For comparison, we put in the numbers for the quartic contribution: 

[(e mi x e^,) • (e m2 x e* m ,\ + (e mi x e m2 ) • (e^, x e^, )] (30) 

Combining the quartic and exchange contributions gives the level pattern (0.3829, —0.1571, —.7658) 
in units of N c g 2 u)o /167T, restoring an anti- ferromagnetic coupling. 

It remains to evaluate the annihilation graph. In this case the currents carry the time depen- 
dence e ziz2luJ ° t , so the propagator is no longer the static Green function. For this graph the currents 
turn out to be 

Jk ~ Attuj 2 2uj r 2 [r X (£mi X Em ^ k Jk ~ & 4ttu 2 2uj t 2 V X ^ X 6 ^ )Jfe ' ^ 
and 

A(x) = JdyJ dte^Du^t-y^Jkti) 
satisfies the inhomogeneous Helmholtz equation 

<-v - .V, = Jl = * («-. * <•;)]<■ <*) 

with cj = 2wo, but it is clearer to write the solution for general u>. Again with the ansatz Ai = 
A(r)[r x (e mi x e m2 )]i, this becomes 

A» + U' + u?A = ^»^, (33) 

with the boundary condition remaining RA'(R) + 2A(R) = 0. Putting these into the formula for 
the contribution to the energy shift gives 

, . r>K 2 r R 

(f\AE\i) annih = N c g 2 5° 5%,— s- / r 2 drA(r)h{ujvr) 2 {e mi x e m2 ) • (e* m , x e* m , ). (34) 

Wq JO 12 

We see immediately that this graph only contributes to the shift of the triplet spin 1 state. This 
is, of course, not surprising because it describes a direct mixing between two gluon and one gluon 
states and the one gluon state is pure spin 1. 

The solution for A with the correct boundary conditions is easily found to be, for general uj, 



.. , 9K 2 uj r R , . 

A \ r ) = ~, 9-7; — / r dr 7i(a;r<J 

v ; 4iruj 2 r2uj J JU <; 



uRy (ujR) - yi{ujR) . 

j/i(wr>) . 7 m JiK 

u>Rj G {u)R) - ji(uR) 



hiuor') 2 , (35) 



where we use the standard notation r< (r>) for the lesser (greater) of the pair r, r'. We want uj = 2ujq, 
but with uj general we should recover the previous case for uj — ► 0, a useful check on numerical 
work. Stripping off the color factors, the annihilation graph contribution to the level shift becomes 



{j\A annih \i} = — — ——2 / zdzji(z) 

16vr y(z Y Jo 



n — \ — n — \^ w z > 



16vr x(zo) 2 Jo 

J o z'dz 'ji(vz')ji{z') 2 {e mi x e m2 ) • (e^, x e^, ). (36) 
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where we have defined rj = uj/ujo — > 2 for this graph. For this value of rj, the numerical evaluation 
of the above expression gives 

</|A ann ^|i) « 0.1430^p(e mi x e m2 ) • (e^, x e^). (37) 
Finally putting all contributions together we arrive at the total level shift matrix 
(m , 1 ,m' 2 \A tota [\m 1 ,m 2 ) R* 

^|^[-0.3829(e mi x e*,) ■ (e m2 x e*,) + 0.2559(e mi x e m2 ) • (e^, x £%)], (38) 

1(37]- 12 12 

which translates to the level pattern (0.3829,0.1289,-0.7658) in units of N c g 2 uJo/16TT for spins 
2, 1, respectively. 

The upshot of our calculations is that the perturbative interactions between nearest neighbors on 
the gluonic chain described by 't Hooft's large N c limit have the correct sign and spin dependence 
to provide a satisfactory pattern of single glueball states. But of course perturbation theory is 
only valid for weak coupling (at short distances), and a near neighbor bond of massless gluons 
requires this coupling to be at least of order 1. As the coupling increases at larger distances, 
perturbation theory will break down, but it is perhaps not too much to hope that the qualitative 
pattern of interactions is not drastically changed. The anti-ferromagnetic character of the spin- 
spin interactions is very reassuring. Given that the bonds do form, even a relatively weak residual 
anti-ferromagnetic interaction should be enough to guarantee that the lowest excitations of a long 
chain have low spin. 

It is noteworthy that this pattern would not hold if the only interaction were due to the A 
term: gluon exchange is essential. Thus, for example, the fishnet diagrams, based on a quartic 



coupling and analyzed in the large coupling light-front studies of [11], cannot contain the essential 
physics of large N c QCD. Those diagrams were singled out by setting up a cutoff model in which 
the strong coupling limit "froze" the P + distribution to be uniform amongst all the gluons in the 
chain. In fact, this frozen P + distribution led to a string picture of the gluonic chain with the 
degrees of freedom of a critical string, whereas QCD must be subcritical (4<26!). To bring into 
a strong coupling expansion the exchange effects necessary for anti-ferromagnetic interactions, a 
more sophisticated cutoff model which allows fluctuations in the P + distribution at large coupling 
must be devised. It is from these fluctuations that we can expect the Liouville world sheet field, 
necessary for subcritical string theory, to emerge. As described in Q the Liouville field should also 
be the "5th" dimension in any alleged QCD/5dString connection. 

Acknowledgements. I would like to thank Joel Rozowsky and Klaus Bering for helpful comments 
on the manuscript. 
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